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Abstract  This  study  systematically  investigated  the  ocean  mixed  layer  responses  to  tropical  cyclone  (TC) 
using  available  Argo  profiles  during  the  period  of  1998-201 1  in  the  northwest  Pacific.  Results  reveal  that  iso¬ 
thermal  layer  (IL)  deepening  and  isothermal  layer  (IL)  cooling  with  evident  rightward  biases  induced  by 
strong  TCs  are  clearer  compared  to  the  weak  TCs.  Likewise,  the  rightward  biases  of  IL  deepening  and  cool¬ 
ing  induced  by  fast  TCs  are  more  obvious  than  that  induced  by  slow  TCs.  The  upwelling  within  TC's  eye  is 
much  stronger  for  the  strong  (slow)  TCs  than  weak  (fast)  TCs.  For  the  strong  and  slow  TCs,  the  TC-induced 
rainfall  reduces  deepening  of  constant  density  layer  (with  its  depth  called  the  mixed  layer  depth,  MLD),  and 
in  turn  increases  the  barrier  layer  thickness  (BLT).  The  initial  BL  prior  to  TC  can  restrict  IL  cooling  more  mark¬ 
edly  under  the  weak  and  fast  TCs  than  under  the  strong  and  slow  TCs.  The  inertial  oscillation  is  stronger 
induced  by  the  strong  (fast)  TCs  than  by  the  weak  (slow)  TCs.  In  addition,  the  most  pronounced  TC-induced 
mixed  layer  deepening  and  IL  cooling  in  July  to  October  climatology  occur  in  the  subtropical  gyre  of  the 
northwest  Pacific  with  enhanced  vertical  diffusivity.  The  maximum  increase  of  isothermal  layer  depth  (ILD) 
and  MLD  is  up  to  5  m,  with  IL  cooling  up  to  0.4°C. 


1.  Introduction 

Oceanic  response  to  tropical  cyclones  (TCs)  has  been  a  hot  topic  due  to  its  importance  for  climate  change, 
ecological  variability,  and  environmental  protection.  Many  studies  show  that  strong  winds  associated  with 
TCs  induce  a  decrease  in  isothermal  layer  temperature  (ILT)  and  a  deepening  of  isothermal  layer  (IL) 
beneath  storms  by  vertical  entrainment  process  due  to  turbulent  mixing.  The  cooling  and  deepening  are 
generally  more  evident  on  the  right  side  of  TC's  tracks  in  the  Northern  Hemisphere  [Price,  1981;  Shay  et  a /., 

1 992;  Chu  et  al.,  1 999,  2000a;  Jacob  et  al.,  2000 ;  Zheng  et  a!.,  2008;  Nam  et  ai,  201 1  ].  Under  storms,  energy 
transferred  from  atmosphere  to  ocean  generates  upwelling,  near-inertial  waves,  and  currents  [Price  et  ai, 

1 994;  Wang  and  Huang,  2004;  Liu  et  ai,  2008;  Jaimes  and  Shay,  201 0;  Chang  et  ai.,  201 2,  2013].  As  TC's  trans¬ 
lation  speed  is  less  than  phase  speed  of  the  first  baroclinic  mode,  oceanic  response  tends  to  be  barotropic, 
geostrophic,  and  cyclonic  gyre  with  upwelling  in  the  storm's  center  [Chang  and  Anthes,  1978;  Chang,  1985; 
Ginis  and  Sutyrin,  1995],  which  weakens  the  IL  deepening  near  the  storm's  centers  [Wu  and  Chen,  2012].  If 
the  TC's  translation  speed  is  faster  than  the  phase  speed  of  the  first  baroclinic  mode,  the  currents  in  the 
wake  become  more  near-inertial  after  the  first  half  inertial  period.  The  upper  ocean  changes,  in  turn,  affect 
the  change  of  typhoon  intensity  through  the  air-sea  feedback  mechanism  [Emanuel,  1991].  The  persistence 
of  sea  surface  cooling  for  1-2  months  after  the  TC  passage  could  potentially  affect  large-scale  atmospheric 
circulation  [Hart  et  ai,  2007].  Several  studies,  based  on  observations  and  modeling,  suggest  that  TC-induced 
energy  input  and  mixing  may  play  an  important  role  in  climate  variability  through  regulating  the  oceanic 
general  circulation  and  its  variability  [e.g.,  Emanuel,  2001 ;  Sriver  and  Huber,  2007;  Liu  et  al.,  2008;  Hu  and 
Meehl,  2009;  Wang  et  al.,  2009;  Fedorov  et  al.,  201 0;  Wang  et  al.,  2014]. 

To  understand  the  interaction  between  ocean  and  TC,  subsurface  measurements  following  TC  are  impor¬ 
tant.  The  Array  for  Real-time  Geostrophic  Oceanography  (Argo)  floats  measure  the  global  ocean  (0-2000  m) 
temperature  and  salinity  even  under  extreme  atmospheric  and  oceanic  conditions  such  as  TC  passage,  and 
thus  provide  a  complete  picture  of  the  upper  ocean  structure  and  variability  in  response  to  TCs  [Roemmich 
and  Owens,  2000].  Statistical  characteristics  of  mixed  layer  responses  to  typhoons  in  the  North  Pacific  are 
obtained  and  examined  by  Park  et  al.  [2005]  using  the  Argo  profiling  float  data  from  2000  to  2003.  The 
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results  indicate  that  the  changes  in  ILT  and  isothermal  layer  depth  (ILD)  are  inversely  correlated,  and  maxi¬ 
mum  cooling  and  deepening  appear  at  midlatitude  of  20°N-30°N,  which  is  mainly  associated  with  strong 
wind,  slow  translation  speed  of  typhoons,  and  shallow  background  ILD.  The  linear  regression  (data  from 
2000  to  2005)  between  the  Argo  profile  and  wind  speed  [Liu  et  al.,  2007]  shows  a  weakly  positive  correlation 
between  the  ILD  change  and  the  wind  speed,  and  implies  that  the  high  wind  speeds  not  only  induce  strong 
vertical  mixing  (increasing  ILD)  but  also  cause  upwelling  (decreasing  ILD).  The  response  of  the  upper  ocean 
to  typhoons  is  reexamined  in  the  northwest  Pacific  Ocean  using  the  Argo  profile  data  (2000-2008)  with  the 
seasonal  variability  removed  [Wu  and  Chen,  2012].  The  result  reveals  a  delayed  response  of  ILT  relative  to 
ILD  and  the  longer  restoring  time  of  ILT.  Above  statistical  analyses  mainly  reveal  the  local  relationship  (one¬ 
dimensional  in  vertical)  between  ILD  and  ILT  from  all  TCs. 

Although  the  two-dimensional  response  to  TCs  in  the  upper  ocean  has  been  studied  [e.g .Jacob  etal,  2000; 
Chu  et  al,  2000a],  statistical  verifications  have  yet  been  conducted  using  the  longer  period  Argo  profile 
data.  Since  the  feedback  of  the  upper  ocean  to  TC  mainly  occurs  during  its  passages,  it  is  necessary  to  study 
the  spatial  structure  of  the  instantaneous  response  of  the  upper  ocean. 

In  regions  of  high  freshwater  input  (i.e.,  heavy  rainfall)  where  significant  salinity  stratification  sets  within  a 
deep  isothermal  layer,  a  barrier  layer  can  appear  between  the  base  of  the  isothermal  layer  and  the  base  of 
the  mixed  layer  (i.e.,  constant  density  layer,  shallower  than  the  isothermal  layer)  to  change  the  thermohaline 
front  [Chu  and  Wang,  2003],  and  to  reduce  the  entrainment  of  the  cool  thermocline  water  into  the  isother¬ 
mal  layer  [e.g.,  Chu  et  al,  2002;  Chu  and  Fan,  201 0,  201 1 ;  Liu  et  al.,  2009].  Usually,  the  TC-induced  sea  surface 
cooling  negatively  feedbacks  to  TC  intensification  [e.g.,  Schade  and  Emanuel,  1999;  done  and  Uhlhorn, 

2003].  The  reduction  of  entrainment  weakens  the  TC-induced  surface  cooling  [e.g.,  Wang  et  al.,  201 1 ;  Bala- 
guru  et  al.,  2012],  which  in  turn  decreases  the  negative  feedback  effect  of  upper  ocean  to  TC.  However,  due 
to  the  observational  limitation,  the  BL  variability  during  TC  has  yet  been  explored  systematically. 

With  increasing  number  of  Argo  profiling  floats  deployed  in  the  world's  oceans  over  the  past  two  decades, 
more  detailed  analysis  on  the  oceanic  response  including  BL  to  TCs  can  be  conducted.  In  this  study,  all  TCs 
during  1998-201 1  are  classified  into  four  types  according  to  their  maximum  sustained  wind  speeds  (weak, 
strong)  and  translation  speeds  (slow,  fast).  The  spatial  and  temporal  variability  of  the  instantaneous  isother¬ 
mal/mixed  layer  responses  to  the  four  types  of  TCs  are  analyzed  and  discussed,  with  the  emphasis  of  the 
barrier  layer  variability  during  TC  passage.  Besides,  the  effect  of  TC  on  the  isothermal/mixed  layer  climatol¬ 
ogy  is  also  investigated  based  on  the  Argo  observation  of  1998-201 1.  The  rest  of  the  paper  is  organized  as 
follows.  Data  source  and  processing  are  described  in  section  2.  Section  3  presents  the  TC  classification  by 
maximum  sustained  wind  speed  or  translation  speed.  Section  4  depicts  the  instantaneous  isothermal/mixed 
layer  responses  to  four  types  of  TCs.  Section  5  discusses  the  temporal-spatial  variability  of  the  isothermal/ 
mixed  layer  responses  to  the  four  TC  types.  Section  6  investigates  the  effects  of  TC  on  the  isothermal/mixed 
layer  climatology.  Summary  and  discussion  are  given  in  section  7. 


2.  Data 


The  6  h  best-track  data  of  TCs  occurring  during  1998-201 1  over  the  northwest  Pacific  are  acquired  from  the 
Joint  Typhoon  Warning  Center,  which  consist  of  position,  maximum  sustained  surface  wind  speed,  and  min¬ 
imum  surface  central  pressure.  Daily  precipitation  with  1°  spatial  resolution  is  obtained  from  Global  Precipi¬ 
tation  Climatology  Project  (GPCP)  [Huffman  et  al.,  2001].  Daily  evaporation  data  with  1°  spatial  resolution 
are  from  OAFlux  [Yu  and  Weller,  2007;  Yu  et  al.,  2008].  The  wind  data  set  is  from  a  new  Cross-Calibrated, 
Multi-Platform  (CCMP)  ocean  surface  wind  product  with  a  spatial  resolution  of  0.25°  X  0.25°  and  a  temporal 
resolution  of  6  h.The  data  set  is  produced  by  combining  all  ocean  surface  wind  speed  observations  from 
SSM/I,  AMSR-E,  and  TMI,  and  all  ocean  surface  wind  vector  observation  from  QuikSCAT  and  SeaWinds.  Fol¬ 
lowing  Oey  et  al.  [2006],  the  wind  stress  is  calculated  from  the  wind  field  using  a  bulk  formula: 


f  0.0012poXl/10|L/10|, 


t  = 


pa  (0 . 00049 + 0 . 000065 \U^\)XU^\UW\, 

pa  (0.001 364+2.34X  1CP5|L/10|  — 2.3X  1CP7|L/10|2)xI/i0|L/io|) 


|U10|  <  llms”1 
11  <  \UW\  <  19ms-1 
19  <  |U10|  <  100ms-1 


(D 


This  formula  is  modified  form  Large  and  Pond  [1981]  to  incorporate  the  limited  drag  coefficient  for  high 
wind  speeds  [Powel  et  al.,  2003]. 
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Figure  1.  Wind  stress  magnitudes  (color  bar)  and  vectors  (black  arrow)  under  weak,  strong,  slow,  and  fast  TCs.  Here  Rmax  is  the  radius  of 
TC's  max  winds  with  47  km  as  the  mean. 


In  this  study,  we  utilize  all  the  available  Argo  profiles,  from  the  Coriolis  Data  Center  (ftp://ftp.ifremer.fr),  for 
1998-201 1  in  the  northwest  Pacific.  In  order  to  keep  data  accuracy,  real-time  and  delayed-mode  quality 
controls  are  performed  using  the  method  of  Wong  et  al.  [2003].  We  choose  the  difference  criteria  to  deter¬ 
mine  ILD  and  mixed  layer  depth  (MLD)  in  the  upper  ocean  from  Argo  profile  data,  which  requires  the  devia¬ 
tion  of  temperature  (density)  from  its  upper  reference  level  value  to  be  smaller  than  a  certain  fixed  value. 
Since  IL  or  ML  is  generally  deeper  than  10  m  in  the  tropical  regions,  the  upper  reference  level  is  chosen  as 
10  m  [Wang  et  al.,  201 1].  The  ILD  is  defined  as  a  depth  at  which  the  temperature  is  0.5°C  lower  than  the 
10  m  depth  temperature  [Chu  et  al.,  2000b,  2002].  The  MLD  is  corresponding  to  the  depth  with  density  dif¬ 
ference  Ao>  from  the  upper  reference  level  (here,  10  m  level),  and  the  A<rf  is  calculated  as  Arrf=0.5^  at 
10  m  depth  [Sprintall  and  Tomczak,  1992].  The  difference  between  ILD  and  MLD  is  defined  as  barrier  layer 
thickness  (BLT). 

3.  TC  Classification 

TCs  are  classified  into  two  groups  either  by  maximum  sustained  wind  speed  or  translation  speed  such  as 
weak  TCs  with  maximum  sustained  wind  speed  <33  m  s-1,  and  strong  TCs  with  maximum  sustained  wind 
speed  >33  m  s-1,  slow  TCs  with  translation  speed  <7  m  s-1,  and  fast  TCs  with  translation  speed  >7  m 
s-1.  Since  the  mean  phase  speed  of  the  first  baroclinic  mode  in  the  northwest  Pacific  is  about  2.88  m  s-1 
during  summer  [Chang  et  al.,  2013],  the  Froude  number  Fr,  which  is  the  ratio  between  the  translation  speed 
of  TCs  and  the  phase  speed  of  first  baroclinic  mode,  is  more  than  2  for  fast  TCs,  and  less  than  2  for  slow  TCs. 

The  wind  stress  vectors  and  magnitudes  calculated  from  CCMP  are  represented  by  the  matched  TC-coordinate 
systems  with  the  unit  vectors  in  the  along-track  and  cross-track  directions.  The  composite  of  the  wind  stress 
vectors  and  magnitudes  under  different  types  of  TCs  are  shown  in  Figure  1.  The  previous  study  of  Me/  et  al. 
[2012]  suggested  that  storm  intensity  correlates  with  translation  speed,  with  hurricanes  of  category  5  moving 
on  averaged  1ms"1  faster  than  tropical  storm.  Figure  1  shows  that  the  wind  stress  of  the  fast  TCs  is  stronger 
than  that  of  the  slow  TCs.  Therefore,  our  result  here  roughly  agrees  with  the  result  of  Mei  et  al.  [201 2].  The  left- 
to-right  asymmetric  distribution  of  the  wind  stress  under  the  fast  TCs  is  clearer  compared  to  the  slow  TCs. 
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4.  Instantaneous 
Response  of  the  Mixed 
Layer 

To  analyze  instantaneous 
response  of  upper  ocean  to  dif¬ 
ferent  types  of  TCs,  a  match-up 
database  is  first  established 
between  the  Argo  profile  pairs 
and  the  TCs  in  the  following 
three  steps:  (1)  For  each  TC,  the 
Argo  profiles  within  500  km 
radius  to  its  center  and  within 
6  h  after  the  passage  are 
selected.  The  time  selected 
increment  of  6  h  equals  the 
temporal  resolution  of  the  TC 
data  to  avoid  a  single  Argo 
profile  matching  to  two  TCs.  (2) 
The  selected  Argo  profile  is 

paired  with  another  profile  taken  before  the  passage  of  the  TC  to  represent  the  variation  induced  by  the  TC. 
The  maximum  distance  between  the  paired  is  restricted  within  50  km  [Wu  and  Chen ,  2012].  This  scale  selec¬ 
tion  is  based  on  a  trade-off  between  the  data  availability  and  the  closeness  of  data  points  to  estimate  the 
local  response.  The  maximum  time  interval  between  them  is  restricted  to  10  days.  (3)  The  selected  Argo 
profile  pairs  and  corresponding  TC  are  matched  up  to  establish  the  first  database.  The  red  dots  in  Figure  2 
show  the  distribution  of  605  Argo  pairs  in  the  database. 


Figure  2.  Selected  Argo  profile  pairs  used  in  section  4  (red  dots)  and  section  5  (blue  dots). 


Argo  profiles  from  the  first  database  in  the  Cartesian  coordinates  also  are  rotated  into  the  matched  TC- 
coordinate  systems.  Figure  3  shows  the  number  of  profile  pairs  within  1 0  times  of  the  radius  of  maximum 
wind  (Rmax)  in  the  cross-track  and  the  along-track  directions.  The  Rmax  with  47  km  is  the  mean  of  the  radius 
of  TCs  maximum  winds  [ Hsu  and  Yana ,  1998].  In  the  most  area,  the  slow  TCs  usually  have  the  most  profile 
pairs  in  all  of  the  four  groups.  In  contrast,  the  fast  TCs  generally  have  the  least  profile  pairs.  The  ILD,  MLD, 
BLT,  and  ILT  can  be  obtained  from  Argo  profiles,  and  their  variability  during  TC  passage  is  calculated  as  the 
difference  between  the  Argo  profile  pairs  to  represent  real-time  responses  of  the  mixed  layer  to  TCs.  Figure 
4  gives  the  statistical  characteristic  of  changes  for  the  ILD,  MLD,  BLT,  and  ILT  under  the  different  types  of 
TCs.  The  uncertainty  is  estimated  by  the  standard  deviation  divided  by  the  root  mean  square  of  the  number 
of  profile  pairs  in  each  bin  [Wu  and  Chen ,  201 2].  The  uncertainties  in  the  inner  core  of  all  TCs  except  for  the 
fast  TCs  are  small.  The  large  uncertainty  away  from  the  TC  center  and  under  fast  TCs  may  be  related  to 
fewer  samples  (Figure  5). 


The  net  freshwater  flux  is  obtained  from  daily  evaporation  and  precipitation  during  the  period  from  the  1 
day  before  TC  passage  to  the  day  of  TC  occurrence,  which  is  shown  by  the  black  contours  in  Figure  4.  The 
rightward  biases  of  the  IL  and  ML  deepening  and  IL  cooling  are  more  obvious  during  strong  TCs  than  dur¬ 
ing  weak  TCs.  The  magnitude  and  extent  of  the  IL  deepening  induced  by  strong  TCs  are  larger  than  those 
induced  by  weak  TCs.  The  maximum  IL  deepening  on  the  right  side  of  the  strong  TC's  center  is  above  1 5  m. 
The  magnitude  of  ML  deepening  induced  by  the  strong  TCs  is  larger  in  the  most  area  except  for  the  area 
with  large  net  freshwater  flux  compared  to  weak  TCs.  The  reduction  of  ILT  on  the  right  side  of  strong  TC's 
center  is  usually  larger  than  that  to  the  weak  TC's  center,  with  extensive  coverage  (Figure  4). 


Similarly,  the  magnitude  of  the  IL  deepening  on  the  right  side  of  the  fast  TCs  is  larger  than  that  on  the  right 
side  of  the  slow  TCs  (Figure  4),  which  may  be  related  to  the  stronger  wind  stress  on  the  right  side  of  the 
fast  TCs  (Figure  1 ).  MLD  changes  are  not  completely  coherent  with  the  ILD  changes.  The  ML  deepening  on 
the  right  side  is  larger  for  the  slow  TCs  than  for  the  fast  TCs.  Like  ILD,  the  stronger  wind  stress  on  the  right 
side  of  the  fast  TCs  supports  that  the  isothermal  layer  cooling  induced  by  the  fast  TCs  is  wider  than  that 
induced  by  the  slow  TCs  (Figure  4). 
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Figure  3.  Number  of  profile  pairs  used  in  section  4  for  the  (a)  weak,  (b)  strong,  (c)  slow,  and  (d)  fast  TCs.  Here  Rmax  is  the  radius  of  TC's  max 
winds  with  47  km  as  the  mean. 


In  the  TC  inner  core  (i.e.,  within  about  100  km  of  the  TC  center),  the  net  freshwater  fluxes  under  the  strong 
and  slow  TCs  are  much  larger  than  those  under  the  weak  and  fast  TCs.  For  strong  and  slow  TCs,  the  corre¬ 
sponding  BLT  in  the  area  with  large  net  freshwater  flux  increases  rapidly.  The  positive  buoyancy  flux  of  the 
TC-induced  rainfall  can  affect  static  stability  of  the  ocean  mixed  layer,  with  the  ability  to  restrain  the  vertical 
mixing.  Thus,  the  TC-induced  rainfall  tends  to  make  ML  shoal  and  in  turn  to  increase  BLT. 

The  impact  of  the  initial  BLT  prior  to  TC  on  the  instantaneous  response  of  ILT  is  analyzed  by  the  linear 
regression  based  on  the  Argo  profile  pairs  within  the  radius  of  200  km  from  TC  centers  (Figure  6).  The  initial 
BLT  is  positively  correlated  with  the  ILT  change  with  the  95%  confidence  level,  indicating  that  the  initial  BL 
can  restrain  the  IL  cooling  [Wang  et  ai,  201 1].  The  correlation  coefficient  between  the  initial  BLT  and  ILT 
change  for  the  strong  TCs  is  smaller  than  that  for  the  weak  TCs.  This  is  because  more  vigorous  vertical  mix¬ 
ing  and  upwelling  induced  by  strong  wind  may  more  easily  overcome  the  haline  stratification  to  entrain 
the  cold  thermocline  water.  Compared  with  the  fast  TCs,  the  longer  residence  of  the  vertical  mixing  and 
upwelling  under  the  slow  TCs  tend  to  support  a  smaller  correlation  coefficient  between  the  initial  BLT  and 
ILT  change. 


5.  Temporal  and  Spatial  Variability  of  the  Mixed  Layer  Response 

To  examine  evolution  of  the  mixed  layer  response  to  TCs,  the  second  match-up  database  is  established 
between  Argo  profile  pairs  and  TCs  by  the  following  step:  (1)  for  each  Argo  profile,  it  will  be  selected  if  there 
is  TC  occurrence  within  the  500  km  away  from  the  profile  and  within  10  day  before  the  observation  of  the 
profile,  and  the  selected  Argo  profile  is  matched  with  the  spatially  closest  TC;  (2)  the  selected  Argo  profile  is 
paired  by  the  similar  method  used  to  establish  the  first  database  in  section  4;  (3)  the  Argo  pairs  and 
matched  TC  are  used  to  establish  the  second  database.  The  blue  dots  in  Figure  1  show  the  distribution  of 
3927  Argo  pairs  in  the  second  database. 

Based  on  the  second  database,  the  spatial-temporal  variability  of  the  mixed  layer  response  under  different 
types  of  TCs  is  examined  with  highly  varying  spatial-temporal  distributions  of  the  profile  pair  number  from 
40  to  5  in  each  bin  (Figure  7).  The  spatial-temporal  distribution  of  ILD,  MLD,  BLT,  ILT  changes  after  TC 
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Figure  4.  Instantaneous  response  of  the  ILD  (m),  MLD  (m),  BLT  (m),  and  ILT  (°C)  during  the  weak,  strong,  slow,  and  fast  TCs.  Black  contour 
indicates  the  averaged  net  freshwater  flux  (mm/d)  from  the  1  day  before  TC  occurrence  to  the  day  of  TC  occurrence.  Here  Rmax  is  the 
radius  of  TCs  max  winds  with  47  km  as  the  mean. 


passage  are  given  in  Figure  8.  For  weak  TCs,  there  are  no  obvious  rightward  biases  of  IL,  ML  deepening  and 
IL  cooling  during  about  1 0  days  after  TC  passage  (Figure  8).  At  the  third  to  fourth  days  after  weak  TC  pas¬ 
sage,  within  the  radius  of  200  km  the  IL  and  ML  deepen,  and  the  IL  cools  (Figure  8).  Such  change  is  more 
evident  under  the  strong  TCs  than  under  weak  TCs.  The  ML  shoaling  within  1 00  km  of  the  strong  TC  centers 
is  much  clearer,  which  may  result  from  the  vigorous  upwelling  in  the  region  of  the  TC  inner  core.  The  associ¬ 
ated  IL  cooling  under  the  strong  TC  has  more  evident  rightward  bias.  The  maximum  IL  cooling  can  reach  up 
to  about  1°C  near  the  strong  TCs  center  (Figure  8). 

Besides,  the  rightward  bias  of  IL  deepening  is  more  evident  induced  by  the  fast  TCs  than  by  the  slow  TCs, 
which  is  attributed  to  the  stronger  left-right  asymmetry  distribution  of  the  wind  stress  under  fast  TCs  (Figure 
1 ).  The  ML  shoaling  near  the  inner  core  of  the  slow  TC  is  clearer  than  that  near  the  inner  core  of  the  fast  TC 
as  the  result  of  strong  upwelling.  Compared  to  the  fast  TCs,  the  temporal-spatial  distribution  of  IL  cooling 
induced  by  the  slow  TCs  presents  a  left-right  symmetry  character  with  the  maximum  cooling  near  the  TCs 
center,  which  indicates  the  important  role  of  upwelling  (Figure  8). 

Previous  studies  showed  that  the  Froude  number  Fr  determines  the  forms  of  the  upper-ocean  response: 
upwelling  or  near-inertial  wave  wakes  [ Geisler ,  1970;  Nilsson ,  1995].  If  Fr>  1,  the  baroclinic  response  with 
the  near-inertial  waves  plays  a  dominant  role,  whereas  ocean  response  tends  to  be  barotropic  response 
with  the  upwelling  in  the  storm  centers  when  Fr  becomes  to  smaller  and  smaller.  Thus,  the  mixed  layer 
response  to  the  slow  TCs  with  Fr  <  2  may  be  dominated  by  the  stronger  upwelling. 
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Figure  5.  Uncertainty  of  ILD  (m),  MLD  (m),  and  ILT  (°C)  changes  during  the  weak,  strong,  slow,  and  fast  TCs,  which  is  calculated  by  the 
standard  deviation  divided  by  the  root  mean  square  of  the  number  of  profile  pairs.  Here  Rmax  is  the  radius  of  TCs  max  winds  with  47  km 
as  the  mean. 

The  uncertainty  is  also  estimated  by  the  standard  deviation  divided  by  the  root  mean  square  the  number  of  pro¬ 
file  pairs  in  each  bin  (Figure  9).  The  uncertainties  of  ILD,  MLD,  ILT  changes  in  most  area  are  small,  with  mean  values 
of  4.0  m,  2.3  m,  and  0.25°C,  respectively.  The  large  uncertainty  is  mainly  associated  with  few  Argo  profile  samples. 

Figure  10  shows  the  radial  dependence  of  the  azimuthally  averaged  ILT  and  ILD  changes  within  10  day  after 
TC  passage  for  the  weak,  strong,  slow,  and  fast  TCs.  ILT  and  ILD  variability  under  the  weak  TCs  is  relatively 
small  (Figure  1 0a).  The  IL  deepening  induced  by  the  strong  TCs  has  a  maximum  of  8  m  at  about  1 50  km  of 
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Figure  6.  Relationships  between  initial  barrier  layer  thickness  (m)  and  ILT  changes  (°C)  within  the  radius  of  200  km  for  the  weak,  strong, 
slow,  and  fast  TC. 
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Figure  7.  Time-space  (radial)  plots  of  the  azimuthally  averaged  profile  pairs  for  the  (a)  weak, 
(b)  strong,  (c)  slow,  and  (d)  fastTCs.  The  positive  (negative)  sign  of  the  horizontal  axis  indicates 
right  (left)  side  of  the  along-track  direction.  The  unit  of  time  is  the  day  after  TCs  passage. 


the  right  side  of  the  along-track 
direction.  It  is  noted  that  the 
strong  IL  deepening  is  also 
found  at  about  150  km  of  the 
left  side.  This  may  be  caused  by 
the  passage  of  TCs  over  the 
ocean  leaving  behind  the 
strong  inertial  oscillations  and 
thus  leading  to  the  formation  of 
convergence  and  divergence 
zones  [ Ginis ,  2002;  Rao  et  a!., 
2010;  Wang  et  ai,  2014].  The 
maximum  IL  cooling  of  about 
0.6°C  induced  by  the  strong 
TCs  is  situated  at  50  km  of  the 
right  of  the  TC  track,  where  the 
IL  deepening  reaches  minimum 
as  a  result  of  the  vigorous 
upwelling  (Figure  10b).  Com¬ 
pared  with  the  ILD  change,  the 
ILT  change  exhibits  a  more 
obvious  asymmetry  with  more 
cooling  on  the  right  side  of  the 
strong  TC  track,  which  may  be 
related  to  the  existence  of  both 
TC-induced  vertical  mixing  and 
upwelling  mechanisms. 


The  ILT  and  ILD  changes  induced  by  the  slow  TCs  show  the  strongly  positive  correlation  with  the  correlation 
coefficient  of  0.77  on  the  95%  confidence  level  (Figure  10c).  Generally,  the  TC-induced  vertical  mixing  can 
make  IL  deepen  whereas  the  TC-induced  upwelling  can  make  IL  shoal.  Besides,  both  mechanisms  can  make 
IL  cool.  Thus,  the  positive  correlation  suggests  a  dominant  role  played  by  the  Ekman  pumping  under  the 
slow  TCs.  However,  the  correlation  between  the  ILT  and  ILD  change  induced  by  the  fast  TCs  is  out-of-phase, 
with  the  correlation  coefficient  of  -0.70  which  is  statistically  significant  at  the  95%  confidence  level.  This 
negative  correlation  under  the  fast  TCs  indicates  the  dominance  role  of  the  mixing  process. 


Time  evolutions  of  the  averaged  ILT  and  ILD  changes  within  the  radius  of  100  km  for  the  weak,  strong,  slow, 
and  fast  TCs  clearly  show  the  inertial  oscillation  induced  by  TC  passage  with  three  minima  of  the  ILT 
changes  within  1 0  day  for  each  TC  type  (Figure  1 1 ).  Comparing  with  the  weak  TCs,  ILD,  and  ILT  changes 
show  larger  amplitude  within  4-1 0  days  after  the  strong  TC  passages,  indicating  that  the  inertial  oscillation 
induced  by  the  strong  TCs  can  persist  for  a  longer  time.  Figure  1 1  b  shows  that  the  IL  cooling  lags  1  day 
behind  the  IL  deepening  after  the  strong  TC  passage.  The  lagged  correlation  coefficient  between  the  ILT 
and  ILD  changes  is  -0.84,  which  is  statistically  significant  at  the  95%  confidence  level.  The  amplitudes  of 
the  ILT  and  ILD  changes  are  smaller  within  1  -5  days  induced  by  the  slow  TCs  than  by  the  fast  TCs,  indicating 
the  stronger  inertial  oscillation  after  the  fast  TC  passages.  The  ILT  and  ILD  sharply  reduce  within  8-1 0  days 
after  the  fast  TC  passages.  Just  as  the  instantaneous  response,  the  relationship  between  the  initial  BLT  and 
the  ILT  change  is  also  analyzed  by  linear  regression  within  10  days  after  TC  passages  (Figure  12).  The  corre¬ 
lation  coefficients  between  the  initial  BLT  and  the  ILT  change  are  all  positive  for  the  weak,  strong,  slow,  and 
fast  TCs.  It  also  proves  that  the  initial  BL  can  weaken  IL  cooling. 


6.  Effect  of  TCs  on  the  Mixed  Layer  Climatology 

Since  TCs  pronouncedly  occur  in  July  to  October,  the  Argo  profiles  are  selected  to  investigate  the  effect  of 
TCs  on  the  mixed  layer  during  July  to  October.  The  Argo  profiles  are  separated  into  two  groups.  The  first 
group  includes  all  the  Argo  profiles  in  July  to  October  during  1998-201 1,  whereas  the  second  group 
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Figure  8.  Time-space  (radial)  plots  of  the  azimuthally  averaged  ILD  (m),  MLD  (m),  BLT  (m),  and  ILT  (°C)  changes  induced  by  the  weak, 
strong,  slow,  and  fast  TCs.  The  positive  (negative)  sign  of  the  horizontal  axis  indicates  right  (left)  side  of  the  along-track  direction.  The  unit 
of  time  is  the  day  after  TCs  passage. 


excludes  the  Argo  profiles  influenced  by  TCs.  The  number  of  Argo  profiles  in  3°  X  3°  grid  cells  for  each  of 
two  groups  is  shown  in  Figures  13a  and  13b,  and  the  difference  between  two  groups  is  given  in  Figure  13c. 
More  than  50  Argo  profiles  are  available  in  each  bin  with  above  1 50  Argo  profiles  near  the  Kuroshio  for 
each  of  two  groups.  The  Argo  profiles  under  the  influence  by  TCs  are  mainly  located  in  the  area  of  1 20°- 
140°  longitudes  and  15°-35°  latitudes  where  TCs  frequently  occur  (Figure  13c). 

Based  on  all  the  Argo  profiles  during  1998-201 1,  the  means  of  ILD,  MLD,  and  ILT  in  3°  X  3°  grid  cells  are  cal¬ 
culated  to  present  the  spatial  distribution  of  climatology  in  July  to  October  (Figure  14).  The  MLD  reduces 
northward  and  southward  from  45  m  in  the  zonal  band  of  10°N-20°N  to  less  than  25  m  at  high  latitudes. 
The  ILT  decreases  gradually  from  low  to  high  latitudes  with  much  higher  than  29°C  at  low  latitudes,  and 
less  than  1 7°C  at  high  latitudes  (north  of  40°N).  The  spatial  distributions  of  climatological  ILD,  MLD,  and  ILT 
obtained  from  the  first  Argo  profile  group  are  consistent  with  the  previous  studies  [e.g.,  Chong  et  al. ,  2009]. 

Figure  15  shows  the  differences  of  climatological  ILD,  MLD,  and  ILT  between  the  first  and  second 
groups.  Figure  16  shows  the  number  of  TC  occurrence  and  the  averaged  maximum  sustained  wind 
speed  in  3°  X  3°  grid  cells  during  1998-2011.  Under  the  TC  influence,  the  IL  and  ML  deepen  up  to 
5  m  in  the  subtropical  gyre  of  the  northwest  Pacific  (125°E-155°E,  20°N-40°N).  The  averaged  maximum 
sustained  wind  speeds  are  much  larger  than  those  in  other  areas,  which  indicate  stronger  TC-induced 
mixing  (Figure  16b).  Actually,  the  most  pronounced  vertical  diffusivity  induced  by  TCs  situates  in  this 
region  [Sriver  and  Huber,  2007].  The  associated  IL  cooling  induced  by  TCs  in  the  TC-active  region  is 
more  than  0.1 2°C,  with  the  maximum  cooling  of  about  0.4°C.  However,  the  IL  and  ML  become  shallow 
and  IL  cooling  becomes  weak,  even  warming  in  the  outer  edge  of  TC-active  region,  for  example,  in  the 
most  area  of  (140°E-165°E,  5°N-25°N)  where  the  probability  of  TC  passage  is  very  low.  The  positive 
wind  stress  curl  in  the  TC-active  region  causes  the  upwelling  and  horizontal  divergence,  and  in  turn 
produce  the  downwelling  and  horizontal  convergence  in  the  outside  edge  of  the  TC-active  region.  Our 
results  also  support  that  the  TC-winds  can  strengthen  climatological  background  Ekman  pumping  [ Jul - 
lien  et  al.,  2012;  Scoccimarro  et  al.,  201 1]. 
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the  weak,  strong,  slow,  and  fast  TCs.  The  positive  (negative)  sign  of  the  horizontal  axis  indicates  right  (left)  side  of  the  along-track  direction. 
The  unit  of  time  is  the  day  after  TCs  passage. 


Distance  (km)  Distance  (km) 

Figure  10.  Radial  dependence  of  the  azimuthally  averaged  ILT  (blue  line,  °C)  and  ILD  (green  line,  m)  changes  within  10  days  after  TC  pas¬ 
sage  under  the  weak,  strong,  slow,  and  fast  TCs.  The  positive  (negative)  sign  of  the  horizontal  axis  indicates  right  (left)  side  of  the  along- 
track  direction.  The  error  bar  with  ±1  standard  deviation  indicates  the  estimated  uncertainty. 
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Figure  1 1 .  Time  evolution  of  the  averaged  ILT  (blue  line,  °C)  and  ILD  (green  line,  m)  changes  within  the  radial  distance  of  1 00  km  under 
the  weak,  strong,  slow,  and  fast  TCs.  The  error  bar  with  ±1  standard  deviation  indicates  the  estimated  uncertainty. 


7.  Summary  and  Discussion 

Based  on  the  historical  data  (1998-2011)  of  Argo  profiles,  TCs,  CCMP  winds,  evaporation,  and  precipitation, 
the  characteristics  of  the  mixed  layer  response  to  different  types  of  TCs  have  been  obtained  in  this  study  for 
the  northwest  Pacific.  The  main  conclusions  are  given  as  follows. 


1 .  The  instantaneous  response  of  the  mixed  layer  to  the  different  types  of  TCs  is  investigated.  Compared  to 
the  weak  TCs,  the  rightward  biases  of  IL  and  ML  deepening  and  IL  cooling  are  more  evident  during  the 
strong  TCs  as  a  result  of  the  more  obvious  rightward  biased  distribution  of  wind  stress  under  the  strong 
TCs.  It  is  demonstrated  from  the  observation  that  the  heavy  rainfall  in  the  inner  core  of  the  strong  and 
slow  TCs  tend  to  make  ML  shoal,  and  therefore  BL  thicken.  Similarly,  the  magnitudes  of  IL  deepening 
and  IL  cooling  on  the  right  side  of  the  fast  TCs  with  stronger  wind  stress  are  larger  than  that  of  the  slow 
TCs.  It  is  more  obvious  that  the  initial  BL  prior  to  TC  passage  restrains  the  IL  cooling  under  the  weak  and 
fast  TCs,  relative  to  the  strong  and  slow  TCs. 
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Figure  12.  Relationship  between  initial  barrier  layer  thickness  (m)  and  ILT  changes  (°C)  under  (a)  the  weak,  (b)  strong,  (c)  slow,  and  (d)  fast 
TC.  The  Argo  profile  pairs  within  the  radial  distance  of  200  km  from  the  TCs  centers  and  10  days  after  the  TC  passage  are  used. 
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Figure  13.  Number  of  (a)  all  the  Argo  profiles,  (b)  the  Argo  profiles  excluding  that  influenced  by  TCs,  and  (c)  the  Argo  profiles  influenced 
by  TCs  in  the  3°  X  3°  grid  cells  during  1 998-201 1 . 


Figure  14.  Climatological  mean  of  (a)  ILD  (m),  (b)  MLD  (m),  and  (c)  ILT  (°C)  from  all  the  Argo  profiles  during  July  to  October  of  1998-201 1 
with  3°  X  3°  resolution. 


2.  The  temporal-spatial  variability  of  the  mixed  layer  response  to  the  different  types  of  TCs  is  explored. 
Comparing  with  the  weak  TCs,  the  variability  of  the  ILD,  MLD,  and  ILT  is  stronger  within  1 0  days  after  the 
strong  TC  passage,  with  strong  upwelling  in  the  inner  cores.  The  rightward  biases  of  IL  and  ML  deepen¬ 
ing  and  IL  cooling  induced  are  more  evident  and  the  upwelling  in  their  inner  cores  induced  is  weaker  by 
fast  TCs  than  by  slow  TCs,  which  may  be  attributed  to  the  stronger  left-right  asymmetry  distribution  of 
wind  stress  and  faster  translation  under  the  fast  TCs.  The  mixed  layer  response  is  dominated  by  the 
Ekman  pumping  for  the  slow  TCs,  whereas  by  the  mixing  process  for  the  fast  TCs.  The  ILT  and  ILD 
changes  within  the  radius  of  100  km  for  all  the  TCs  show  the  inertial  oscillation  with  larger  amplitudes 
under  strong/fast  TCs  than  under  weak/slow  TCs.  Furthermore,  the  IL  cooling  lags  1  day  behind  the  IL 
deepening  within  the  radius  of  100  km  from  strong  TC  centers. 

3.  The  TC  effect  on  mixed  layer  climatology  during  July  to  October  is  studied.  The  TCs  act  to  deepen  the  cli¬ 
matological  ILD  and  MLD,  and  cool  the  ILT  in  the  area  with  frequent  TC  passages.  The  most  pronounced 
IL,  ML  deepening,  and  IL  cooling  occur  in  the  subtropical  gyre  of  the  northwest  Pacific  (125°E-155°E, 
20°N-40°N)  where  strong  vertical  diffusivity  is  induced  by  TCs.  The  maximum  IL  and  ML  deepening  is  up 
to  5  m,  with  IL  cooling  up  to  0.4.  Nevertheless,  it  is  also  found  that  the  IL  and  ML  become  shallow  and 
the  IL  cooling  becomes  weak,  even  warming  in  the  outer  edge  of  TC-active  region.  It  is  mainly  caused  by 


Figure  15.  Differences  of  (a)  the  ILD  (m),  (b)  MLD,  and  (c)  ILT  (°C)  between  all  the  Argo  profiles  and  the  Argo  profiles  under  no  influence  of 
TCs  during  July  to  October  of  1 998-201 1  with  3°  X  3°  resolution. 
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Figure  1 6.  Horizontal  distributions  of  (a)  total  number  of  6  h  occurrence  of  TCs,  and  (b)  averaged  maximum  sustained  wind  speed  (m/s)  in 
3°  X  3°  grid  cells  during  1998-2011. 


Latitude  (°) 


Figure  17.  Latitudinal  distribution  of  averaged  (a)  MLD  and  (b)  ILT  changes  (m)  from  Argo  profile  pairs  under  no  influence  of  TCs  (red)  and 
under  the  influence  of  TCs  (blue). 


the  upwelling  and  horizontal  divergence  due  to  positive  wind  stress  curl  in  the  TC-active  region,  and  in 
turn,  to  balance  the  water  mass,  the  downwelling  and  horizontal  convergence  occur  in  the  outside  edge 
of  the  TC-active  region. 
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When  we  analyze  the  mixed  layer  changes  between  after  the  TC  passage  and  before  the  TC  passage, 
the  statistical  method  may  have  some  uncertainty,  which  mainly  result  from  sampling  errors  and  oce¬ 
anic  background  variability.  The  sampling  error  has  been  given  and  discussed.  Here  the  uncertainty 
caused  by  the  background  variability  is  given  and  discussed.  All  Argo  profiles  of  1998-2011  are  divided 
into  two  groups  according  to  whether  they  are  influenced  by  TCs  or  not.  For  each  group,  Argo  profiles 
are  also  paired  by  the  similar  method  in  section  4  and  the  changes  of  the  MLD  and  ILT  are  calculated 
based  on  Argo  profile  pairs.  The  averaged  MLD  and  ILT  changes  under  non-TC  and  TC  influence  are 
given,  respectively  (Figure  17).  Though  the  oceanic  background  MLD  changes  at  10°N-30°N  and  ILT 
changes  at  30°N-50°N  are  larger  than  those  at  other  latitudes,  they  are  far  less  than  the  TC-induced 
changes.  Thus,  the  results  in  the  study  can  represent  the  characteristics  of  the  mixed  layer  response  to 
the  TCs  with  a  permissible  uncertainty. 
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